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1 . INTRODUCTION 


The Space and Missile Systems Organization of the Department of Defense 
(DOD) has recently issued a Request for Proposal (RFP) for the design of an 
Interim Upper Stage ( I US ) to be used in conjunction with the Space Shuttle 
for transporting STS payloads to orbits beyond the Shuttle performance 
envoi ope. 

The IUS design is to employ expendable solid rocket motors (SRM's) for 
major translational maneuvers. Although the question is not addressed spe- 
cifically in the RFP, it is expected that the IUS design will not provide 
for guidance-commanded SRM thrust termination such as might conceivably be 
effected by chemical quenching or thrust chamber venting. A summary of prob- 
able IUS design characteristics, extracted from the aforementioned RFP, is 
included in the Appendix. 

The main purpose of this report is to discuss techniques of trajectory 
design, maneuver execution, and stage loading that are compatible with the 
use of SRM's which, once ignited, must burn to propellant depletion. It is 
anticipated that some Shuttle payloads will use non-IUS solid propellant kick 
stages; therefore this' subject is also pertinent to Shuttle flights other 
than those involving the use of the IUS. 

The SRM utilization techniques can be divided into two major categories: 
(1) those in which the stage performance is adjusted to match the a V require- 
ments of a preselected trajectory, and (2) those in which the trajectory is- 
designed to match the aV capability of the stage(s). These two categories’ 
are discussed separately in Sections 2 and 3. 
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2. ADJUSTMENT OF SRM STAGE PERFORMANCE 
TO MATCH- TRAJECTORY AV REQUIREMENTS 


Within this category of techniques it is assumed that a nominal trajectory 
has been selected to deliver a particular payload to a specific orbit and that, 
for each major velocity impulse, an SRM has been identified which is capable 
of delivering more than the required A V. The problem then is one of reducing 
the effective aV delivered by the SRM stage so as to match the trajectory 
requirement. 

These techniques allow the use of conventional trajectory design methods 
and software; i.e., those in which velocity impulse magnitudes are not con- 
strained to be equal to specific values. Aside from any question of compu- 
tational difficulty, it is not desirable to impose equality constraints on 
the velocity impulses during trajectory design because doing so removes, de- 
grees of freedom which may be needed to satisfy environmental or operational 
constraints. 

2.1 STAGE BALLASTING 

For a given specific impulse I sp , propellant loading W p , and burnout 
gross weight Wg^, the ideal velocity increment that will be delivered by 
the SRM is uniquely determined by the rocket equation 

4V ideal * 9 ! SP ^ ' (') 

Ballasting increases the burnout weight Wg^, thereby reducing the aV capability 
of the stage. Although useful in some instances, this technique is generally 
undesirable because it involves a double weight penalty. The weight of the 
excess propellant and the weight of the ballast necessary to counteract it 
must both be transported into orbit by the Shuttle. 


The velocity increment that would be realized by thrusting in a fixed 
direction in gravity-free space. 
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2.2 PROPELLANT OFFLOADING 

The aV capability of a given SRM stage can also be reduced by offloading 
propellant, thereby reducing Wp in Equation (1). Weightwise, this is more 
efficient than ballasting. The magnitude of the weight advantage over bal- 
lasting will depend on whether or not it is necessary to replace the offloaded 
propellant with inert material so as to maintain the structural integrity of 
the SRM, which in turn depends mainly on the propellant grain configuration. 

A practical problem associated with SRM propellant offloading is that 
it must be accomplished at the manufacturing facility, well in advance of the 
flight date. Therefore this technique is intolerant of changes in payload 
weight or flight objective after the SRM has been assembled. 

2.3 POWERED-FLIGHT ATTITUDE MANEUVERS 

Neither of the two preceding techniques affords, within themselves, the 
capability to compensate for trajectory dispersions detected during flight. 

If the upper stage guidance and control system is capable of controlled at- 
titude maneuvers during the SRM burn, then such compensations can be made by 
varying the orientation of the thrust vector. Given sufficient maneuver- 
ability, gross reductions of the effective SRM aV can be attained. The basic 
idea is to waste the excess aV produced by the SRM by thrusting in a "non- 
optimum" direction. Many variations of this technique are possible; two of 
the simpler ones are discussed in Sections 2.3.1 and 2.3.2. 

2.3.1 Dog-Leg Steering 

Figure 1 illustrates the geometry, in inertial velocity space, of the 
dog-leg maneuver. This is essentially an "external aV" maneuver type, in. 
which the objective is to achieve a resultant a V vector (as sensed by onboard 
accelerometers) supplied as a target value by premaneuver computations. 
Operationally, the SRM thrust line would be aligned with the inertial direc- 
tion of Leg 1 prior to ignition. After ignition, the thrust direction would 
be held fixed until the proper time (as determined by accelerometer output) 
for initiating a constant-rate turn that would align the thrust vector 
with Leg 2. After achieving proper alignment on Leg 2 (on a bearing that 
passes through the coordinates of the desired resultant aV), the thrust 
direction would again be held inertially fixed until SRM burnout occurred. 


3 




For a given SRM aV capability greater than the required resultant aV, 
there are obviously many different combinations of turning rate and Leg 1/ 

Leg 2 orientation that would nominally achieve the desired result. The 
turning rate will of course be limited by permissable load factors and by 
the available directional control authority and control accuracy. Assume 
for the moment that a rate has been determined. This still leaves open the 
question of how the direction of Leg 1 or of Leg 2 should be chosen. If the 
orientation of either leg is specified, the orientation of the other is 
uniquely determined by the SRM AV capability. 

Since the total aV delivered by the SRM (consequently, the length of 
Leg 2) will always be uncertain to some extent, it is believed that the best 
policy for dog-leg steering would be to align Leg 2 with the "least critical 
direction" in velocity space. By definition, this is the direction in which 
a velocity variation of a given magnitude will have the smallest possible 
undesirable effect on the ensuing trajectory. The definition of the least 
critical direction would depend on the maneuver objective. For instance, 
if the maneuver were calculated to intercept a moving target, this might be 
the direction (easily calculable from a state transition matrix or from solu- 
tion of Lambert's problem) in which a small velocity variation would not cause 
the target to be missed, but instead would only change the time of intercept. 
In general, the least critical direction can be defined by the vector cross 
product of the gradients, in velocity space, of the two most critical 
postmaneuver trajectory parameters. 

2.3.2 Circular-Arc Steering 

The geometry of the circular-arc steering algorithm (Reference 1) is 
illustrated in Figure 2. Unlike the dog-leg method, its use is not restricted 
to "external aV" maneuver types. It can also be used in conjunction with 
guidance logic in which the required velocity ^ is periodically updated 
during the maneuver (e.g., by Lambert solutions based on the navigation sys- 
tem's estimate of the current state vector). At any instant, the desired 
angle a between the thrust vector and the velocity- to- be gained vector tg is 
calculated from the equation 

IfJ 

sm a _ 1 G 1 
a ~ V ’ 
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where V is the guidance system's estimate of the remaining A V capability 
of the SRM. As illustrated in Figure 2, V and |$ G | correspond to the arc 
length and chord lengths respectively, of a segment of a circle. A potential 
advantage of this steering policy is that it tends to minimize the turning 
rate, although in practice it may be necessary to freeze the thrust direction 
near the end of the burn to preclude excessive rates that could result from 
underestimating V. 

It should be noted that the thrust direction is not completely defined 
by Equation (2). The steering angle a can be measured in any plane that 
contains 1L. And for any permissable plane, either of two opposite angular 
directions may be chosen. To minimize the deleterious effects of SRM AV 
uncertainties, probably the best policy would be to choose the steering plane 
and angular direction that would align the burnout thrust vector as closely 
as possible with the "least-critical" direction discussed in Section 2.3.1. 
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3. DESIGN OF TRAJECTORIES TO MATCH SRM 
STAGE AV CAPABILITIES 

For SRM stages incapable of powered-flight attitude maneuvers, it may 
be necessary to design trajectories with velocity impulse magnitudes con- 
strained to be equal to specific stage AV capabilities. Even if ballasting 
and/or propellant offloading techniques are used to adjust SRM stage per- 
formance to the requirements of a preselected nominal trajectory, such a 
trajectory design capability may be needed during flight operations so that 
maneuvers can be retargeted to compensate for dispersions. 

In Sections 3.1 through 3.4, various classes of upper-stage maneuver 
targeting and trajectory design problems are discussed. Emphasis is placed 
on the nature and types of the solutions that are possible when AV magnitudes 
are constrained to specific values, as deduced from analytical considerations. 

A conic-impulsive trajectory model, wherein coasting flight segments are ap- 7 
proximated by two-body motion and powered flight segments are approximated 
by instantaneous velocity increments, is used in all cases for analytical 
purposes. Such a model is adequate in most cases for conceptual design and 
for feasibility assessment. Analytic solutions based on such modeling are 
also valuable for high-fidelity precision maneuver targeting. They can serve 
either to generate "first guesses" for conventional iteration processes that 
use state transition matrices, or to provide ‘both the first guesses and updated 
independent variable values (based on trial trajectory integrations) in a 
recursive iteration loop similar to that described in References 2 and 3. 

3.1 SINGLE-IMPULSE ORBIT MODIFICATION 

The problem types to be considered in this section are those where it 
is desired to attain, with a single velocity increment of fixed magnitude, 
a postmaneuver orbit having specific characteristics. In the general case 
where the impulse can be applied at any point on the premaneuver orbit, it 
may be possible to satisfy as many as three mutually independent postmaneuver 
orbit constraints. However, the discussion in this section will be limited 
to the restricted case where the maneuver position and the premaneuver veloc- 
ity are fixed; e.g., by specifying either the time of the impulse or the 
angular position of the impulse point in the premaneuver orbit. With this 
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type of restriction there are only two degrees of freedom (represented, 
say, by two angles that define the direction of the velocity increment), 
and therefore no more than two mutually independent postmaneuver orbit 
constraints can be satisfied simultaneously. 

For analysis purposes, consider an inertial Cartesian coordinate system 
in which the Z-axis is oriented in a direction opposite to that of the posi- 
tion vector of the impulse point, the Y-axis is oriented in a direction op- 
posite to that of the angular momentum vector, and the X-axis completes a 
right hand orthogonal triad. The impulse position and premaneuver vectors 
can then be represented by 


and 

A ^ A 

v_ = i X + k Z (4) 

As illustrated in Figure 3, the locus of all. possible postmaneuver velocity 
coordinates is represented by the surface of an "SRM capability sphere 11 in 
velocity space. The equation of the capability sphere is 

(* - X Q ) 2 + Y 2 + (z - 1 Q ) 2 = V 2 , (5) 

where V represents the magnitude of the velocity increment produced by 
the SRM, and where the postmaneuver velocity vector is represented by 

v = iX + jY + kZ . (6) 

In general, the specification of one postmaneuver orbit constraint 
defines a surface in velocity space. The specification of two constraints 
defines one or more curves in velocity 'space representing the intersection(s) 
of two surfaces, provided of course that it is physically possible to satisfy 
the two constraints (i.e., that the constraint surfaces do in fact intersect 
each other). The maneuver targeting problems discussed in Sections 3.1.1 
through 3.1.5 involve geometric solutions to find the points, if they exist, 
at which such curves intersect the surface of the SRM. capability sphere. 


9 





The problem types to be discussed are representative of those commonly 
encountered, but by no means do they encompass all possible single-impulse 
targeting problems. 

3.1.1 Orbit Circularization 

Requiring the postmaneuver orbit to be circular is the same as requiring 
the eccentricity e to be zero. In general, assignment of a specific value 
to e defines a surface in velocity space just as if it were any other scalar 
orbit parameter. At e = 0, however, the constraint surface degenerates into 
a curve. Specifically, the locus of velocity coordinates that satisfy this 
constraint is a circle whose equations are 

Z=0 (7) 

and 

X 2 + f = u/r o , (8) 

where y is the gravitational parameter of the central body. Substituting 
Equations (7) and (8) into Equation (5) yields the partial solution 

(•' /r o + *o 2+2 o 2 - v2 )/< 2 > i o> (9) 

"k 

provided X Q f 0 . The remainder of the solution is found by substituting 
the numerical value of X from Equation (9) into Equation (8) and solving 
for the two possible values of ?, thus: 

Y* = ±\/(u/r 0 ) - $ 2 , (10) 


< 

If = 0 (which is unlikely, since it could only occur if the premaneuver 

velocity vector were either zero or oriented in a vertical direction), then a 

solution exists only if y/r +X 2 +Z 2 ~V 2 =0. 

J 0 0 0 
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where the subscript i takes on the value 1 when the minus sign is used 
and the value 2 when the plus sign is used. If the radicand is negative 
in Equation (10), then of course there is no physical solution to the 
targeting problem. 

In general, the roots of Equation (10) will not be zero. This means 
than an orbital plane change is generally necessary if the orbit is to be 
circularized with a velocity increment of fixed magnitude. The magnitude 
and direction of the plane change is defined by the wedge angle W, where 

sin W = V V + f 2 (11) 


and 


cos U = X/ v X + Y 


4 


( 12 ) 


The wedge angle is positive when the maneuver produces a positive rotation 
of the orbit plane about the Z-axis. 

3.1.2 Achievement of Specified Apogee and Perigee Altitudes 

Specifying the altitude of either apsis (apogee or perigee) of the 
postmaneuver orbit is a special case of a more general type of constraint 
which requires that the flight path angle p (the angle between the position 
vector and the velocity vector) have a specific value at a given radial 
distance r-j- from the center of attraction. Since the general form is appli- 
cable to other targeting problems (e.g., calculating deorbit maneuvers) , the 
general equation of the constraint surface will be developed before invoking 
the special conditions that apply to apsis altitude constraints. 

For the trivial case where r x = r , the constraint surface is that of 

T o 

a cone whose equation is 

Z 2 sin 2 e T = {% 2 + 'Y' 2 ) cos 2 3j . (13) 
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For the non-trivia! case (i.e., where ry f r Q ), the conservation of energy 
and angular momentum requires that 

X 2 + t 2 + 2 2 - 2vi/r Q = Vy 2 - 2p/r T (14) 

and that 

r Q 2 (X 2 + ? 2 ) = r T 2 v T 2 sin 2 By , (15) 

where Vy is the velocity magnitude in the postmaneuver orbit at distance ry 
from the center of attraction. Substitution of the expression for Vy 2 from 
Equation (14) into Equation (15) yields, after some algebraic manipulation, 
the equation 


(X 2 + Y 2 )[l - r Q 2 /(r T sin By) 2 ] 



+ 



(16) 


If r T < r , then Equation (16) defines a hyperboloid of one sheet whose 
intersection with the X-Z plane is shown in Figure 4. If ry > r , then the 
geometric figure represented by Equation (16) depends on sin By. When 
sin By < (r Q /ry) the figure is a hyperboloid of two sheets that degenerates 
into two parallel planes at sin By = (>” 0 /ry). When sin By > (r Q /ry), the 
figure is an ellipsoid. The intersections of the constraint surfaces with 
the X-Z plane for these latter three cases in shown in Figure 5. In all 
cases, the surfaces are symmetrical about the origin and about the Z-axis. 

The applicable surface in three-dimensional velocity space for each case can 
be visualized as that which is swept out by the appropriate curve in Figure 4 
or Figure 5, as it is rotated about the Z-axis. 

Now consider the special case where ry represents the desired distance 
of one of the apsides (either perigee or apogee) from the center of attrac- 
tion. If the subscript ct is used to designate apsidal conditions in lieu 
of the general subscript T, then r T is replaced by r^, cos By is replaced by 
0, and sin By is replaced by 1 in the preceding equations, since the flight 
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Figure 5. Intersection of r T , 3y Constraint Surfaces 
with X-Z Plane (r T > r ) 





path angle g must be equal to 90 degrees at either apsis. Again, taking 
up first the trivial case where r^ = r , Equation (13) becomes 

l 2 = 0 , (17) 

which is of course the equation of a horizontal plane in velocity space. 

For the non-trivial case (where r^ f r Q ), Equation (16) can be simplified 
to the form 


(k 2 + Y 2 )(l + r o /r a ) _ ^2 

+ - V r a T 


(18) 


Note that the first term in Equation (18) is always postive, and that the 
second term is positive or negative depending on whether r^ is greater 
than or less than r . Physically, of course, if r is greater than r n it 
must represent the apogee distance r^, and if it is less than r Q it must 
represent the perigee distance r p . For the general case where r p < r Q < r A , 
the two. equations 


(X 2 + Y 2 )(l + r o /r A ) ^ 2 2 

(2 Wry + (2y/r o )(l - rjr A ) " 1 


and 


(X 2 + ? 2 )(r o /r p + 1) ^ 2 

(Wyi (2u/r Q )(r o /r p - l) ' 1 


(19) 


(20) 


must be satisfied simultaneously. Equations (19) and (20), respectively, 
define an ellipsoid and a hyperboloid of one sheet whose intersections with 
the X - Z plane are illustrated in Figure 6. The intersections of the 
ellipsoid and hyperboloid with each other form two parallel circles in 
velocity space whose equations are 


Z 


a 


+ 



(1 - r o /r A )(r o /r P ‘ 
(r 0 /i-p + r 0 /r fl ) 


(21) 
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PERIGEE 

CONSTRAINT 


Figure 6 


Z 

Intersection of Apogee and Perigee Constraint 
Surfaces with X-Z Plane (r < r < r A ) 




and 


X^ + Y = 


(2y/r 0 ) 

(r 0 /rp + r o /r A ) 


( 22 ) 


and where the subscript & is assigned a value of 1 when the minus sign 
is used and a value of 2 when the plus sign is used. The two circles can 
be visualized as those which are swept out' by the points of intersection be- 
tween the ellipse and one branch of the hyperbola shown in Figure 6, as they 
are rotated about the Z-axis. 

Given r , r^, and r p , Equation (21) can be used to compute the two 
numerical values of Z that will satisfy the targeting constraints. These 
values can be substituted into Equation (5), along with the expression 


2 _ (2y/r Q ) 2 

" Cr o /r p + W ’ 


(23) 


from Equation (22), to obtain two linear equations in X of the form 


*o 2 + - V 2 + < 2 >^0>/< | -c/ r P + Va 1 - t . 

2 *o 


(24.) 


where i takes on the values 1 and 2. Finally, the targeting solution 
is completed by substituting the numerical values of X-j and X 2 back into 
Equation (23), obtaining 



(2y/r o ) 

-v (r 0 /r p + r 0 /r A ) * 


(25) 


Again, the subscript a takes on the values 1 and 2, and for each value 
of 4 the subscript m is assigned a value of 1 when the minus sign is 
used in Equation (25), and a value of 2 when the plus sign is used. Of 
course, only real values of m are acceptable. 
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If the radicand in Equation (25) is positive for both values of £, 
then there are four physical solutions of the maneuver targeting problem. 

As in the case of the orbit circularization problem, the values of m 
will not be zero in general , which means that an orbital plane change gen- 
erally will be necessary if specific values of apogee and perigee altitude 
are to be attained with a single velocity increment of fixed magnitude. 

3.1.3 Achievement of Specified Orbit Period and Apsis Altitude 

Requiring the postmaneuver orbit to have a specified period P defines 
a sphere in velocity space whose equation is 

f +. Y 2 + Z 2 = (2y/r o ) - (2wu/P) 2/3 . (26)* 


If it is desired also to achieve a specified apsidal distance r , then 
Equations (18) and (26) must be satisfied simultaneously. The loci of 
solutions (if any exist) are two parallel circles whose equations are 


1 


£ 


+ 



(r Jr 
O a 


- 1 ) 





(r <A</ - 1 

. (r 0 /r a )2 


(27) 


and 


• ? •? 

r + r = 


(r Jr )(2 u /r ) - (2WP) 2/3 


(r Jr )' 
o a 


(28) 


where the subscript is assigned a value of 1 when the minus sign is 
used and a value of 2 when the plus sign is used. 

Given r , r , and P, Equation (27) can be used to compute the two 
numerical values of Z that will satisfy the targeting constraints. These 
values can be substituted into Equation (5), along with the expression 


?2 ' (y.HjWrp) - (2WP) 2/3 _ 

< r A / 


( 29 ) 


*If (2y/r ) - (2 ttJP) 2;/3 < 0, of course there is no physical solution 
to the problem. 
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from the Equation (28) to obtain two linear equations in X of the form 


* 0 2 + (Z st ' iof_ + [(yr a )(2|i/r 0 ) - (2 »u/P) 2/3 ] /(yr tt ) 2 - V 2 (30) 
— JY . 


where & takes on the values 1 and 2. The targeting solution is 
completed by substituting the numerical values of X-j and X 2 back into 
Equation (29) to obtain 



( r n/Tj( 2 V>/0 - (2-rrp/P) 2 / 3 


(r_/r )‘ 
' o a' 



(31) 


A's before, the subscript i takes on the values 1 and 2, and for each 
value of a the subscript m is assigned a value of 1 when the minus 
sign is used in Equation (31) and a value of 2 when the plus sign is used. 
And again, of course, only real values of m represent acceptable solutions. 

3.1.4 Achievement of Specified Orbit Period and Wedge Angle 

Imposing an equality constraint on the wedge angle W (the angle 
through which the orbit plane is rotated by the velocity impulse) defines a 
half-plane in velocity space whose equation is 

X sin W = Y cos W , (32) 

where X is required to have the sign of cos W and ? is required to have the 
sign of sin W. The intersection of this half-plane with the sphere defined 
by Equation (26) is a semicircle whose radius is equal to 

■J (2y/r 0 ) - (2irp/P) 2/3 . 

First consider the case where cos W = 0. Equation (32) then yields 

X- 2 = 0 , (33) 
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which can be substituted into Equations (5) and (26) to obtain 


f * a- i 0 ) z • v 2 



(34) 


and 

? 2 + Z 2 = (2y/r Q ) - (Ztfu/P) 273 , (35) 

where of course the right hand sides of both equations must be greater than 
or equal to zero if a real solution is to exist. Expanding Equation (34) 
and then subtracting it from Equation (35) yields the linear equation 

■ * 0 2 + Z 2 - V 2 + (2p/r ) - (2irp/P) 2/3 

l = _o o . . o (36) 

2 *o 

which can be used to obtain a numerical value of Z. This numerical value 
can be used in the equation 

Y = Jv 2 - X Q 2 - (Z - 1 Q ) Z sin W (37) 

to complete the solution for the special case of cos W = 0. Equation (37) 
is derived from Equation (34) and the condition imposed on Equation (32) 
regarding the sign of 

Now consider the case where cos W f 0. Equation (32) yields 

? = X tan W , (38) 

which can be substituted into Equations (5) and (26) to obtain 

(31 - X Q ) 2 + X 2 tan 2 U + (1 - Z Q ) 2 = V 2 (39) 



Z Q = 0, then a solution is possible only if 
+ Z Q 2 - V 2 + (2p/r Q ) - (2 ttp/P) 2/3 = 0. 
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and 


X 2 + X 2 tan 2 W + Z 2 = (2y/r Q ) - (2 ttp/P) 2/3 . (40) 

Expanding Equation (39) and subtracting it from Equation (40) yields 

' Q 1 - X 0 X = 1 Q t , (41) 

where 

2Q, = (2p/r 0 ) - (2pp/P) Z/3 - V 2 + X Q 2 + Z Q 2 . (42) 

If Z Q = 0, then Equation (41) can be used to obtain a numerical value 
for X, thus 

X = Q-,/X 0 . (43) 

If the computed value of X does not have the same sign as cos W, then 
there is no solution to the problem. If the numerical value of X does 
have the same sign as cos W, then it can be substituted into Equation (38) 
to obtain In this case* the two permissable values of 1 are found by 
substituting the numerical value of X into Equation (40), yielding 

2y/r 0 ) - (2TTP/P) 2 / 3 - (X/cos W) 2 , (44) 

where a is assigned the value 1 when the minus sign is used in Equation 
(44) and the value 2 when the plus sign is used. 

If Z f o, then the solution is found by squaring both sides of Equation 
o 2 

(41) and substituting an expression for Z from Equation (40) into the result, 
thus: 

q/ - 2Q 1 X 0 X + £ 0 2 X 2 = Z o 2 [ (2y/r 0 ) - (2 ttu/P) 2/3 

- X 2 (1 + tan 2 W) ] . (45) 
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Equation (45) can be arranged in the quadratic form 

A ] X Z + A 2 X + A 3 = 0 , 


which has the roots 


where 



4A 1 A 3 



A t = X Q 2 + (Z o /cos W) 2 , 
A 2 - - 2 Ql * 0 , 

and 


A 3 = " 1 o L (2p/r o } " ^^/ p ) 2/3 J . 


( 46 ) 


(47) 


(48) 

(49) 


(50) 


and where a is assigned the value 1 when the minus sign is used in 
Equation (47) and the value 2 when the plus sign is used. Only those real 
values of X £ that have the same sign as cos W satisfy the wedge angle con- 
straint. The acceptable values of X £ (if any) can be substituted into 
Equations (38) and (41) to obtain the associated values of ? and Z. 

3,1.5 Achievement of Specified Apsis Altitude and Wedge Angle 

This targeting problem requires the simultaneous satisfaction of 
Equations (5), (18), and (32). The conditions associated with Equation (32) 
regarding the sign of X and of f, as delineated in Section 3.1.4, must 
also be satisfied. 

Again, consider first the case where cos W = 0. Substitution of 
X 2 - 0 from Equation (33) into Equation (18) yields 

? 2 [l - (r 0 / r a )2 ] + = (^o) O - V r o ) • (51) 
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Equation (34), which resulted from substituting ft = 0 into Equation (5), 
can be subtracted from Equation (51) to obtain 

2Z n ft = (, rjr ) 2 ft 2 - 2Q ? , (52) 

o o a 2 

where 


2Q 2 = V 2 - ( X Q 2 + ft/ j - (2 P /r o )(l - r Q /r a ) . (53) 

If ft = 0, Equation (52) and the sign conditions associated with Equation 
(32) yield immediately 


ft = [/2Q 2 j (r Q /r a )J sin W . (54) 

There is no solution to the problem if Q 2 <_ 0. Assuming Q 2 > 0, the numerical 
value obtained from Equation (54) can be substituted into Equation (51) to 
solve for ft, thus: 

\ -- ±v (1 - V r «> [ ( 2 >‘/ r o) - U + V.) ] (55) 

where a is assigned a value of 1 when the minus sign is used and a value 
of 2 when the plus sign is used. 

If ft f 0, the solution is obtained by squaring both sides of Equation 

0 # o 

(52). Then an expression for ft, obtained from Equation (51), can be 
substituted into the result to produce 

, B ] ft 4 + B 2 ft 2 + B 3 = 0 (56) 


which has the roots 




B 2 " 4B 1 B 3 


/2B 1 


(57) 
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where 


B 1 = (r 0 / r a )4/4 ’ (58) 

B 2^o 2 [’ - W 2 ] - WH • ' 59 > 

and 

B 3 ■ 0 2 2 - i 0 2 (2w/r 0 )(l -r 0 /r a ) , (60) 

and where £ is assigned the value 1 when the minus sign is used' in 
Equation (57) and the value 2 when the plus sign is used. For each of 
the real positive values of t (if any) obtained from Equation (57), the 

out-of-plane postmaneuver velocity component is given by 

• ( 6i > 

The radial velocity components are found by substituting the computed values 
of % into Equation (52). 

Now, consider the case where cos W f 0. Substitution of Equation (38) 
into Equation (18) yields 

X 2 (1 + tan 2 W)(l + r 0 /rj j2 

(ZSTFJ (2 M /r 0 >n - TJFJ - ’ (62) 

from which Equation (39) can be subtracted to obtain 

Z Q Z = Q 3 % Z - X Q X - Q 2 , (63) 

where 

2Q 3 = (r 0 /r a ) 2 (l + tan 2 W) (64) 

and where Q 2 is given by Equation (53). 
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If Z Q = 0, then Equation (63) reduces to 

Q 3 X 2 - * 0 X - Q 2 =. 0 (65) 

which has the roots 

V (*o±VV + 4Q 2 (! 3) /(2 ' 3 3> • ( 66 > 

where % is assigned the value 1 or 2 depending on whether the minus 
sign or the plus sign is used in Equation (66). Only those roots of Equation 
(65) that are real and that have the same sign as cbs W are acceptable. For 
the acceptable values of X^, the corresponding values of are obtained from 
Equation (38). The postmaneuver radial velocity components are obtained from 
Equation (62), thus: 

*l.n, ’ ±\l^/r 0 W - r 0 /r a ) - [l - [rjrf ] (tycos H) 2 (67) 

where, for each of the two possible values of i , the subscript m is assigned 
a value of 1 or 2 depending on whether the minus sign or the plus sign is 
used in Equation (67). 

If 1 f 0, both sides of Equation (63) are squared and an expression 
for Z is obtained from Equation (62) and substituted into the result. This 
yields 


C ] X 4 + C 2 X 3 + C 3 + C 4 X + C 5 = 0 , (65) 

where 

C 1 = Q 3 2 , (66) 

C 2 = " 2Q 3 , (67) 

C 3 = * 0 2 - 2Q 2 Q 3 + (2 q /cos W) 2 [1 - (r 0 /r a ) 2 ] , (68) 

C 4 = 2Q 2 * 0 , (69) 
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and 


c 5 ■ v" - * 0 Z ( 2 ‘ I / 1 'o)< 1 - V r «> • (70) 

The numerical values X (where z = 1, 2, 3, 4) which satisfy Equation (65) 
can be computed by using an explicit algorithm which exists for finding the 
roots of a quartic equation. Of the four possible roots, only those which 
are real and which have the same sign as cos W are acceptable. Given any 
acceptable value X £ , the corresponding value Y £ can be obtained from Equation 
(38) and 2 can be obtained from Equation (63). 

3.2 TWO- IMPULSE ORBITAL TRANSFER WITHOUT PHASING CONSTRAINT 

The orbit of any satellite can be defined completely by six scalar orbit 

ic 

constants, such as the following classical set : 

p semilatus rectum 

e eccentricity 

i inclination 

n right ascension of the ascending node 

u argument of perigee 

tp time of perigee passage. 

That is to say, given appropriate numerical values for these six scalars, 
the position and velocity of the satellite can be calculated at any arbitrary 
time t. 

The first five quantities listed in the preceding paragraph define the 
size and shape of the orbit and its orientation in an equatorial inertial 
reference frame. They provide sufficient information to calculate satellite 
position and velocity vectors as functions of a geometric variable such as 
the true anomaly f. The sixth quantity (t p ) is a phasing constant whose 
value must be known to calculate the satellite position and velocity as a 
function of time. 


This set is completely definitive for all orbits except those involving 
rectilinear motion, where p = 0 and e = 1 and consequently the orbital 
energy cannot be calculated. In such a case (which is not considered to 
be of interest here), the value of some other parameter such as the 
semimajor axis a must be known. 
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The problem to be considered in this section is that of transferring a 
spacecraft from a fully defined initial orbit to a target orbit whose size, 
shape, and orientation (p, e, i, ft, to) are given as input but whose phasing 
constant (t p ) is not specified. Furthermore, the transfer is to be effected 
by two velocity increments, at least one of which is required to have a 
specific magnitude. 

If the positions of both impulse points and the orientation of the 
transfer orbit angular momentum vector are known, then there exists a "singly- 
constrained 11 analytic solution to determine explicitly just which transfer 
trajectories (if any) will satisfy a aV constraint at one end or the other 
of the transfer arc. Often it will be necessary to satisfy arbitrary con- 
straints on velocity increment magnitude at both ends of the transfer arc. 

No explicit solution is known for the latter problem. However, doubly- 
constrained solutions can be found by a straightforward procedure which, 
involves "scanning" singly-constrained solution families. The scanning is 
accomplished by systematically varying the values of two independent variables 
that jointly define the two impulse locations and the direction of the trans- 
fer orbit angular .momentum vector. Depending on the overall nature of the 
trajectory design problem, it may be preferable when generating the singly- 
constrained analytic solutions to satisfy the AV constraint on either the 
final or the initial impulse. Assume for the sake of argument that it is 
chosen to satisfy the constraint on the initial impluse. Then it is a simple 
matter, each time a singly-constrained solution is found, to calculate the 
required final AV as an output variable. Since every member of the solution 
family satisfies the a V constraint on the initial impulse, the solutions to 
the doubly-constrained problem can be found (if they exist) by examining the 
output data to determine which combinations of independent variable values 
produce the desired final a V. 

Of course there are many other transfer trajectory parameters besides 

the magnitude of the unconstrained velocity increment which can be computed 

\ 

as output from the singly-constrained solution, and which may have to be 
considered in the selection of a transfer trajectory. Given an array of 
data resulting from a systematic scan of two independent variables, general- 
purpose software is available for generating a display which depicts the 
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behavior of any selected set of output variables over the entire two- 
dimensional scan space (Reference 4, pp. 12-13). The ability to utilize 
such a capability to solve this type of trajectory design problem is highly 
desirable. The key to successful utilization lies in the identification of 
a pair of scalar independent variables Y and 7 such that a systematic 
variation of their values will generate all possible solutions of interest. 

3.2.1 Selection of Scan-Control Variables 

In general, the planes of the initial orbit and the target orbit will 
be inclined to one another at some angle a. As illustrated in Figure 7, 
the position of the initial impulse can be defined by the angle A j , which 
is measured from a relative node (say, the ascending node of the initial 
orbit on the target orbit plane) in the direction of orbital motion in the 
initial orbit. The position of the final impulse can be defined by the angle 
Ay, measured from the same relative node in the direction of motion in the 
target orbit plane. The transfer angle e is measured from the initial 
impulse point to the final impulse point, in the direction of motion in the 
transfer orbit. The transfer orbit makes the wedge angles Wj and lip respec- 
tively, with the initial orbit plane and the target orbit plane. The wedge 
angles are measured from the appropriate (initial or target) orbit plane to 
the transfer orbit plane, and are positive in the direction of clockwise 
rotation when viewed from a vantage point above the appropriate impulse 
point, looking downward toward the center of gravitational attraction. 

The positions of the impulse points obviously can he defined by specifying 
values for Aj and Ay. However, this combination of variables would not be a 
good choice as the only method of scan control because their systematic vari- 
ation would not generate the solution families which are associated with 
nodal transfer (i.e., where Aj = 0° and Ay = 180° or Aj = 180° and Ay ='360°). 
In a nodal transfer case the number of possible transfer orbit planes is 
infinite. Consequently it is necessary to specify the value of some addi- 
tional input quantity, such as the wedge angle at one of the impulse points, 
in order to generate a unique solution. 

There is no ambiguity at any point regarding the orientation of the 
transfer orbit plane if, say, Aj and Wj are chosen as scan control variables. 
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Figure 7. Geometry of Orbital Transfer Problem 


Again, this would not be a good choice as the only method of scan control 
because (for instance) when Xj = 0 and = a, the number of possible loca- 

tions for the second impulse is infinite and some other input quantity such 
as Xj must be input to generate a unique solution. 

It turns out, unfortunately, that there is no single pair of independent 
variables whose systematic variation will generate all of the solutions which 
might be of interest in this type of trajectory design problem. However, all 
possible solutions can be generated by separately scanning at most two sets 
of independent variables, such as those designated as Set No. 1, 2, 3, and 
4 in Table 1.. As indicated in the table. Sets 1 and 2 apply when the singly- 
constrained analytic solution is to satisfy an equality constraint’ on the 
magnitude of the initial impulse, and Sets 3 and 4 apply when an equality 
constraint is to be applied to the final impulse magnitude in the analytic 
solution. 

For any one set of independent variables, there are certain combinations 
of values which will produce "manifold ambiguties" (i.e., situations in which 
there are infinite number of possible solutions) of the type previously dis- 
cussed, wherein a whole family of solutions is obscured by the particular 
choice of independent variables. The conditions which characterize these 
situations are shown in Table 1, along with suggested "resolutions" (i.e., 
arbitrary conditions to be imposed for the purpose of generating one unique 
solution from the obscured family). The true resolution of any manifold 
ambiguity requires a separate scan of the appropriate alternative variable 
set, which allows the obscured family to be surveyed in its entirety. 

In addition to the manifold ambiguities associated with particular 
combinations of independent variable values, there is a "twofold ambiguity" 
associated with every combination of values. In the cases where Aj and Ay 
are specified, this arises from the possibility of two opposite orientations 
of the transfer orbit angular momentum vector. The suggested resolution in 
such cases is to require the cosine of the wedge angle at the constrained 
impulse point to have a particular sign, as defined by an auxiliary input 
variable H that is required to have the value + 1.0 or - 1.0 throughout 
any particular scanning operation. In the cases where Aj and Wj or where 
Ay and Wy are specified, the twofold ambiguity arises from the fact that 
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Table 1. Scan Control Data 


INDEPENDENT 

VARIABLE 

SET 

NO. 

(ISET) 


IMPULSE 
AT WHICH 
AV CONSTRAINT 
IS APPLIED 


INITIAL 


INITIAL 


DEFINITION OF 
SCAN CONTROL VARIABLES 


RESOLUTION 
OF TWOFOLD 
AMBIGUITY 


MANIFOLD AMBIGUITY 

CONDITION 
FOR OCCURRENCE 

RESOLUTION 

1 sin W T = 0 

sin 2 e =0 

sin e =0 

2 2 
sin Wj = sin a 

sin Wj = 0 

sin 2 o=0 

sin 0 =0 

2 2 
sin W T = sin A 















the transfer angle e may have either of two values which differ by 180°. 

The suggested resolution in such cases is to have the aforementioned variable 
H define the sign of cos e. 

The equations shown in Table 1 for the resolution of twofold ambiguities 
are designed such that, in each case, the subset of transfer orbits associated 
with a value of H = - 1.0 will generally require lower aV magnitudes than the 
subset associated with a value of H = + 1 .0. This increases the probability 
that all possible solutions for a particular problem can be generated by using 
only one of the two allowable values of H. 

3.2.2 Selection of Reference Node 

In the general case where sin a » 0, it makes little difference which 
of the relative modes between the initial and target orbit planes is used as 
a reference for measuring Xj and Xy. Either the ascending or the descending 
node of one orbit on the plane of the other might be chosen, so long as it 
is used consistently. Hov/ever, when sin a « 1 (i.e., when the orbits are 
nearly coplanar), special considerations may have to be taken into account 
if discontinuities are to be avoided in the output data between separate 
scans over neighboring orbits, or between successive discrete solutions in 
an iterative precision targeting process. 

When the initial orbit and the target orbit are coplanar, the relative 
node line is of course undefined. The orbits can be coplanar only if 
ij = iy or if ij = 180° - iy. In all cases but the very special one where 
sin ij = sin iy = 0, it is also necessary that = &y when ij = iy or that 
= £2y +180° when ij = 180° - iy. Because gravitational harmonics generally 
produce different values for the nodal regression rates sir and ,Qy , perfect 
coplanarity between satellite orbits is at most a fleeting circumstance in 
the real world of orbital flight operations. 

Figure 8 illustrates a typical case of the type under consideration 

where, for simplicity, the plane of the target orbit is assumed to be Station- 
er 

ary. If perfect coplanarity occurs at the time t , then between the times 
t - e and t + e there is a discontinuous change of 180° in the inertial 
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Figure 8. Node Shift Geometry for Nearly Coplanar Orbits, 









locations of both the ascending and the descending relative nodes. If either 
of these nodes were used as the reference for measuring Xj and Xy , then 
clearly there would also be, within the same time interval, discontinuous 
changes of 180° in the inertial directions of impulse position vectors 
defined by fixed values of Xj and Xy. 

It can be seen by inspecting Figure 8 that the northerly and southerly 
(as distinguished from the ascending and descending) relative nodes approach 
the positions defined by u = 90° and u = 270°, respectively, as e (hence 
also sin a) approaches zero, where u is the argument of latitude in either 
the initial or the target orbit. It follows that the discontinuities just 
described can be removed if Xj and Xy are referenced to, say, the northerly 
relative node when sin A f 0 and to the position defined by u = 90° when 
sin a = 0. 

Although the convention described in the preceding paragraph is to be 
preferred when the initial and final orbits are in the neighborhood of 
coplanarity, unfortunately it can not be used generally for all orbit com- 
binations. To do so could produce discontinuities, similar to those pre- 
viously described, whenever the relative nodes happen to lie in the vicinity 
of the equator. Since the choice of convention is necessary to avoid dis- 
continuities between rather than within separate scans or discrete solutions, 
there is no basis for making such a choice within the logic of the analytic 
solution itself. Therefore it is necessary to provide an input flag IREF 
which can be given either of two values by the trajectory analyst, specifying 
whether the reference for measuring Xj and Xy is to be the northerly relative 
node or the ascending relative node. 

3.2.3 Calculation of Position and Velocity Vectors in the Initial and 
the Target Orbits 

The purpose in this section is to define equations for calculating 
position and velocity vectors in the initial and the target orbits as func- 
tions of the independent variable set identification flag ISET and the as- 
sociated scan control variables and _Y which are defined in Table 1, the 
auxiliary input variable H, and the reference convention flag IREF. These 
computations are necessary to generate inputs to the transfer orbit velocity 
vector equations which will be defined in Section 3.2.4. Only the case where 
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the aV constraint is to be applied to the initial impulse (ISET = 1 or 2) 
will be addressed specifically. The logic and equations for the other case 
(ISET = 3 or 4) can be obtained by replacing ISET with ISET + 2, changing 
the sign of the transfer angle e, and interchanging the subscripts I and 
T or the words "initial" and "target" wherever they appear in the following 
discussion. 

To facilitate subsequent computations, it is convenient to define the 
geocentric inertial components of a set of unit vectors relating to the 
initial and target orbits. In general for any orbit the unit vector pointing 
in the direction of its ascending node on the equator is given by 


Q = 


cos £2 
sin si 


0 '• 


The unit vector 


J = 


- sin Si sin i 
cos Q sin i 

- cos i 


(74) 


(75) 


points in a direction opposite to that of the angular momentum, and the unit 
vector 


U = Si x 5 . (76) 

is directed toward the point of maximum northerly declination in the orbit. 
The unit vectors 


K = - si cos us - U sin to 


(77) 


and 


I = J x K (78) 

A A A 

complete an orthogonal triad (I, J, K) which defines the local vertical 
coordinate system at nerigee just as the triad (i, j, k) used in Section 3.1 
defines the local vertical coordinate system at an arbitrary point in the orbit. 
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The cosine of the relative inclination between the two orbits is given 


by 


cos A - J-j* • Jj 


(79) 


and the sine by 


sin a = v 1 - cos^ A , (80) 

where Jj and Jy are obtained by substituting the appropriate -elements of the 
initial and target orbits, respectively, into Equation (75). If sin A = 0, 
then Xj and Xy will be measured from a reference defined by the unit vector 


N = Uj 


(81) 


If sin a f 0, then the unit vector 




(82) 


defines the direction of the initial orbit's' ascending node on the target 
orbit plane. This is the reference for measuring Xj and Xy when IREF = 1. 
When IREF = 0, then the "nearly coplanar" convention is to be used and it 
is necessary to test the sign of the third component of the vector defined 
by Equation (82). If this component is negative, then Equation (82) has 
defined the southerly relative node, and the sign of all three components 

A 

of N must be reversed (thereby defining the northerly node) before 
proceeding with the rest of the calculations. 

Now, it is necessary to compute the triad of unit vectors ^ i Q , j Q , k Q j 
which define the local vertical coordinate system at the constrained impulse 
point. These are given by the equations 

k Q = ( 5j x N ) sin X - N cos X , (83) 



( 84 ) 
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and. 


1* = Jr X k 

0 10 


( 85 ) 


where, of course, the scan control variable Y represents Aj for IREF = 1 or 2, 
The sine and cosine of the true anomaly in the initial orbit at the constrained 
impulse point can then be obtained from the equations 


sin f = k x K t 
ool 


( 86 ) 


and 


cos f = k • K t 
ool 


(87) 


The position and velocity vectors at the constrained impulse point are then 
given by 



k rt r 
o o 


( 88 ) 


and 






(e-j- + cos f Q ) + Kj 


sinf 0 ] , 


where 


r o = Pl^ 1 + e I cos f 0 } ‘ 


The local vertical components of v Q are given by the equations 



and 


1 


o 




(89) 


(90) 


(91) 


(92) 
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When ISET = 1, then of course it follows from Table 1 that 


sin Wj = sin Y 


( 93 ) 


and 


cos Wj = cos X ■ 


The equation 


J x = j 0 cos Wj - i 0 sin Wj 


' (94) 


(95) 


defines a unit vector which points in a direction opposite to that of the 
transfer orbit's angular momentum vector. The cosine of the wedge angle at 
the unconstrained impulse point then can be obtained from the equation 


cos Wj = Jj • . (96) 

A trial value for the sine of this angle (whose sign may have to be reversed, 
depending on the results of subsequent tests) is then computed by use of the 
equation 


sin W T = Y 1 - cos 2 W T . (97) 

If sin Wj = 0, then there is a manifold ambiguity regarding the location 
of the unconstrained impulse, which is resolved arbitrarily by setting 

kp = - H k Q , (98) 

where H is the auxiliary input variable (described in Section 3.2.2) which 
must have the value + 1.0 or - 1.0. The unit vector kp points in a direction 
opposite to that of the position vector at the "free" (i.e., unconstrained) 
impulse point. 
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If Equation (97) yields a value of sin Wy / 0, then a trial value of 
kp is computed from the equation 

kp = ^ 5-j. x 3^ j /sin W-p . (99) 

If the quantity H. ^ k Q ♦ kp j is negative, then the signs of all three 
components of kp and the sign of sin Wy must be reversed before continuing 
the calculations. 

Having fixed the values of sin Wy and kp, the remaining two unit vectors 
of the local vertical triad at the unconstrained impulse point are obtained 
from 


and 




x kp 


( 100 ) 


( 101 ) 


and the sine and cosine of the true anomaly in the target orbit are obtained 
from the equations 


sin fp = kp x Ky • Jy 


( 102 ) 


and 


cos f = kp • Ky 


(103) 


The position and velocity vectors in the target orbit are given by 





k F r F 


(104) 


and 




+ COS fp) + Ky 


sin fp 



(105) 
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where 


r p = P T /(1 + e T cos fp) . (106) 

The local vertical components of-Vp are obtained from the equations 

X F = v p • i F ' (107) 

and 

^ = v F • k F . (108) 


Finally for the case where ISET = 1, the sine and cosine of the transfer 
angle are given by 


and 


A A A 

sin 6 = k F x k Q • J x (109) 


cos e = kp ■ k Q . ’ (110) 

When ISET = 2, then Y_ = \.p and the value of kp is given unambiguously 
by the equation 

kp = ( J T x N ) sin 7 - N cos T . (Ill) 

The values of jp, ip, sin fp, cos fp, rp, Vp, rp, Xp, and ip can then be 
obtained from Equations (100) through (108), and the value of cos 0 from 
Equation (110). A trial value of sin 0 is computed by use of the equation 


sin 0 W 1 - cos^ 0 . (112) 
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If sin 6=0, there is a manifold ambiguity regarding the direction of the 
transfer orbit's angular momentum vector, which is resolved arbitrarily by 
setting 


J x = - H J I 


(113) 


If sin 0^0, then a trial value of is computed by using the equation 


= ( k F x K) /Sin 


(114) 


If the quantity H (jj • ) is greater than zero, then the signs of sin e 

and all three components of must be reversed before continuing the 


calculations. 


After having fixed the values of sin e and J^, the sines and cosines 
of the wedge angles are computed from the equations 


sin = x J)( * , (115) 

A A 

cos Wj = Jj • , (116) 

sin Wy = J-j- x 0j( • kp , (117) 

and 

cos Wy = J T • J x . (118) 

3.2.4 Calculation of the Transfer Orbit Velocity Vector at the Constrained 
Impulse Point 

As in the last section, only the case where the aV constraint is applied 
to the initial impulse will be addressed specifically here. The logic and 
equations for the other case (where the aV constraint is applied to the final 
impulse) can be obtained by changing the sign of e and interchanging 
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the subscripts I and T, the words "initial" and "target," and the words 
"premaneuver" and "postmaneuver," wherever they appear in the following 
discussion. 

For the purposes of calculating the transfer orbit, it is convenient 
to rotate the local vertical coordinates at the constrained impulse point 
through the wedge angle Wj as illustrated in Figure 9, by means of the 
transformation equations 

T = X cos Wj + Y sin Wj , (119) 

S = - X sin Wj + Y cos Wj , (120) 

and 

R = - Z . (121) 

The R and T axes lie in the transfer orbit plane, and the S axis is 
perpendicular to it. The inverse transformation is defined by the equations 

X = - S sin Wj + T cos Wj , (122) 

Y = S cos Wj + T sin Wj , (123) 

and 

Z = - R . (124) 

In the rotated coordinate system, the postmaneuver velocity coordinates 
are designated R, S, and f, and the equation of the SRM capability sphere 
takes the form 


* • O t • 9 • • . 9 ? 

(R - R q ) 2 + (S - S o ) 2 + (T - T q ) 2 = V 
which is equivalent to Equation (5) with 



(125) 


(126) 
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Figure 9. 


Rotation of Local Vertical Coordinates into 
Transfer Orbit Plane 
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S 0 = - X 0 sin Wj , (127) 

and 

t 0 = X Q cos Wj . (128) 

2 • 2 

If V < S , then the capability sphere does not intersect the required 

0 Z * 2 

transfer orbit plane, and no solution is possible. If V > S q , then the 

capability sphere does intersect the required transfer orbit plane to form 

a circle in velocity space which is defined by the equations 

S = 0 (129) 

(r-R 0 ) 2 + (T-f 0 ) 2 = * 2 -S 0 2 . (130) 

The existence of a solution depends on whether or not this circle intersects 
a hyperbola, also lying in the transfer orbit plane, v/hose equation is to 
be developed shortly. 

Figure 10 is a combined position and velocity vector diagram pertaining 
to the problem at hand. The equations for calculating the transfer angle e 
and the two position vectors r Q and r p were defined in Section 3.2.3. The 
chord vector c can be expressed as 

c = c c , (131) 

where 

c = | r p - f Q | (132) 

and 

c = ( r p - r Q ) /c . (133) 
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The sine and cosine of the angle 2<(>, which will be used later, can then be 
obtained from the equations 


sin 2<j> = k Q x c • 


(133) 


and 


cos 2<f> = k • c , (134) 

o 

where, again, the equations for k Q and 5^ were defined in Section 3.2.3. 

Battin has shown (Reference 5, pp. 104-107) that in general the locus 
of coordinates for all possible transfer orbit velocity vectors (i.e., those 
velocity vectors which, if supplied to a spacecraft at position r , would 
cause its orbit to pass through the position r-p) is a hyperbola whose asymp- 
totes coincide with the chord vector c and with the extension of the position 
vector 'r . The equation of this hyperbola can be developed from the equation 
for the semilatus rectum of the transfer orbit 


and the polar equation 


PX 


P X = (t r//u 


f r F (1 - cos e) 


(135) 


(136) 


o T r - T r c cos e + R r p sin e 


from Reference 6. Equations (135) and (136) can be combined and rearranged 
to yield 

T 2 ( r Q 2 - r Q r F cos e | + f R r Q r F sin e = y r F (1 - cos e) . (137) 


The law of sines and the law of cosines, applied to the triangle shown in 
Figure 10, yields 



and 


2 2 2 
- r r r cos e = I c - r - r r 


o F 


6 = ( 


)/2 , 


which can be substituted into Equation (137) to obtain 


r 


[( 


c + '•o - V 


) /2 ] + T 


T R c r Q sin 2<j) = y rp ( 1 - cos e) 


From the law of cosines we also obtain 


2 2 2 \ 

c + r 0 - fp j /2 = c r Q cos 24 > 


which produces 

"H cos 24 + f R sin 24 = (y/r Q ) (rp/c) (1 - cos e) 


when substituted into Equation (140). 

Now we rotate the R - T coordinate axes through the angle 4 > as 
illustrated in Figure 11 by use of the transformation equations 

p = R cos 4 > - T sin <j> , 

and 


x = R sin 4 > + T cos 4 > . 

The inverse transformation equations 

R = p cos 4 > + x sin 4 > , 


and 


T = - p sin 4 > + t cos 4 > , 

when applied to Equation (142), produce 


.2 

T 


COS 




4 > = (u/r 0 )(r F /c)(l - cos e) 


(139) 

(140) 

(141) 

(142) 

(143) 

(144) 

(145) 

(146) 

(147) 
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after some algebraic manipulation and the substitution of appropriate 
trigonometric identities. 

The right hand side of Equation (147) is always greater than zero 
except for the very special case where cos e = 1. Such a case is of little 
practical interest since it involves vertical rectilinear motion in the 
transfer orbit. Therefore, for our purposes we will assume that no solution 
exists when cos e = 1 . When cos e f 1, Equation (147) can be written in 
the form 


•2 2 
T cos j>_ 


•2 . 2 

, p si n «j> _ -I 

Xp/r 0 ) ( r p / c ) ( 1 - cos- ej " (y/r 0 )(r F /c) d - cos "ej " 
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which can- be recognized as the equation of a hyperbola whose center lies at 
the origin of coordinates and whose major axis coincides with the -r-axis. 

Now we return to Equation (130), or rather to its equivalent'. in the 
p - t coordinate system, which is 

(p - ; o ) 2 + U - i 0 ) 2 - V 2 - S o 2 (149) 

where 

p o = R 0 cos + “ T 0 sin * (150) 

and 

t q = R 0 sin <p + T 0 cos <f> . (151) 

The solutions to the singly-constrained orbital transfer problem (if any 
such exist) must satisfy Equations (147) and (149) simultaneously. 

For the special case where sin <j> = 0 (which occurs when e = 180°), 
Equation. (147) yields immediately 

t = (y/r 0 )(r F /c)(l - cos e) . (152) 
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There is no ambiguity regarding the sign of x. It must be positive; 
otherwise the direction of the transfer orbit's angular momentum vector 
(which, at this point, has already been fixed by the equations in Section 
3.2.3) would be reversed. The numerical value obtained from Equation (152) 
can be substituted into Equation (149) to obtain the two solutions 



where a takes on the value 1 when the minus sign is used in Equation (153) 
and the value 2 when the plus sign is used. 

A value of cos $ = 0 can occur only if cos 0=1, in which case we have 
already decided to assume that no solution exists for the targeting problem. 
For the more general case where neither sin <j> nor cos <j> is equal to zero, 
the expression 


t 2 = p 2 tan 2 <j> + (u/r 0 )(rp/c) (1 - cos 9)(1 + tan 2 <j>.) (154) 

can be obtained from Equation (147) and substituted into Equation (149) to 
produce the equation 

2 t q x = p 2 (1 + tan 2 <ji) - 2p Q p + (155) 

where 

Q 4 = Po 2 + S 0 2 + “ y2 + ^/ r 0 Kr F / c ) 0 - cos e) (1 + tan 2 *) . (156) 

If 4 = 0, then Equation (155) reduces to the quadratic form 

p 2 (1 + tan 2 <j>) - 2p Q p + = 0 (157) 


If the radicand is negative in Equation (153), then of course there is no 
real solution to the targeting problem. 
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which has the roots 


(158) 

where % is assigned the value 1 when the minus sign is used in Equation 
(158) and the value 2 when the plus sign is used. If the roots of Equation 
(157) are real, then they can be substituted into Equation (154) to obtain 

= \j(> 9 Z tan 2 $ + (y/r Q )(r F /c) ( 1 - cos e)(l + tan 2 <j>) (159) 

where, again, r must be positive to maintain the proper direction of the 

. X* 

angular momentum vector. 

If Tq f 0,'then both sides of Equation (155) are squared, and Equation 
(154) is again substituted into the result to obtain the quartic equation 


p 4 + Dg p 2 + D 4 p + Dg = 0 , (160) 

where 

D-j = (1 + tan 2 <f»} 2 , (161) 

°2 = " 4 ^o ^ + tan2 ^ » O 62 ) 

D 3 = 4p Q 2 + 2Q 4 (1 + tan 2 <j>) - 4 t q 2 tan 2 * , (163) 

D 4 = - 4Q 4 p Q , (164) 

and 

°5 = Q 4 2 " 4 ^o 2 ^/ r 0 )(r F /c )(l - cos 9)(1 + tan2 +) * {165) 


p 2 - Q 4 0 + tan 2 <f>) 


(1 + tan^ <j>) 
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Of course, there may be 0, 2, or 4 real roots of Equation (160). Their 
values (if they exist) can be found by using an explicit algorithm for solving 
a quartic equation. If the real roots are designated p , then the corres- 
ponding values t can be obtained from Equation (155), thus: 

X/ 


T *> = 


p 0 2 (1 + tan 2 <j>) - 2p o p £ + _Q 4 


2t 


(166) 


Again, negative values for x £ are unacceptable and must be rejected along 
with the root that produced the negative value when substituted into 
Equation (166). 

3.2.5 Designation of the Desired Solution and Calculation of the 
Transfer Orbit Parameters 

In the preceding section we have seen that there may be as many as four 
distinct real solutions characterized by transfer orbit velocity coordinate 
pairs (p , t ) at the constrained impulse point, where the subscript i can 
have the value 1, 2, 3, or 4. It has already been pointed out that t must 

X/ 

be positive In order to preserve the proper direction of the angular momentum 
vector. In addition, it is necessary that p & + < (2y/r Q ) whenever 

ISET = 1 or 2 and p > 0 or whenever ISET = 3 or 4 and p < 0. This latter 

•V Xt 

restriction is necessary to insure that the transfer orbit, extended in the 
appropriate direction of motion, will actually pass through the position 
defined by rp. 

To obtain a unique solution from the logic described in Section 3.2.4, 
it is necessary for the trajectory designer to assign an- integer value 
between 1 and 4 to an input flag M to designate the particular solu- 
tion that he desires. The value of N must remain constant during any 
particular scanning operation or iterative targeting process, and the selec- 
tion procedure within the analytic solution logic must be such that the output 
data generated by a fixed value of N will vary continuously (at least as 
nearly as possible) as functions of the scan control variables J and T. 

This probably can best be accomplished by ordering all of the acceptable 

coordinate pairs (p n , T n ) such that p n _-j <_ p n < when ISET = 1 or 2 and 

such that p -i > p > p when ISET = 3 or 4. This has the effect of 
n- 1 — n — n+1 
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ordering the acceptable solutions such that At n _-j £ At n <_At n+ .j, where At 
is the flight time which must be spent in the transfer orbit. If there is 
an Nth solution in the ordered set, then it should be used to generate the 
output data. Otherwise, a "no solution" flag should be set to indicate the 
non-existence of the Nth solution, and further computations should be bypassed. 

If the Nth solution does exist, then the local vertical velocity 
components (R, f) in the associated transfer orbit at the constrained impulse 
point can be calculated by substituting p .j and into Equations (145) and 
(146). The transfer orbit's semilatus rectum p x can then be obtained from 
Equation (135), and the eccentricity e x can be computed from the equation 

e x [(p x /r 0 ) - l] 2 + R 2 (p x /m) • (167) 

The sine and cosine of the true anomaly in the transfer orbit at the- 
constrained impulse point can be obtained from the equations 

sin f oX = ^ a/ /e X ( 168 ) 

and 

cos f QX = C(p x /r 0 ) -1] /e x , (169) 

provided e x f 0. If = 0, f x can be assigned an arbitrary value of zero. 

The "local vertical" unit vectors at perigee of the transfer orbit are 
defined by the equations 

h ’ i k o ) cos f ox'- s1n f 0 x (170 > 

and 

*x = (h x O s1n f ox + K cos f oX • < 171 > 
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The sine and cosine of the true anomaly in the transfer orbit at the 
unconstrained maneuver point can be obtained from the equations 

sin fp^ = sin f ^ cos e + cos -f ^ sin e (172) 


and 


cos f 


FX 


cos f ^ cos 


0 + sin f ^ sin e 


(173) 


where the upper sign is to be used when ISET = 1 or 2 and the lower sign 
when ISET = 3 or 4. The velocity vector in the transfer orbit at the 
unconstrained impulse point then can be computed by use of the equation 

*FX w/p X [*X (e X + cos f FX> + s1n f FX ] 1 (174> 

which can be used in the equation 



(175) . 


to compute the magnitude of the unconstrained velocity increment. 

3.3 TWO-IMPULSE ORBITAL TRANSFER WITH PHASING CONSTRAINT 

The problem type to be considered in this section is sometimes referred 
to as the "two-impulse orbital rendezvous problem," since of course to effect 
a rendezvous it is necessary for the active vehicle to adjust all six of its 
orbit constants (p, e, i, a, m, tp) to match those of the target, spacecraft. 

In the true rendezvous problem, where the target is an actual physical object, 
usually there are a number of very stringent constraints relating to natural 
lighting conditions at the time of intercept and to the direction and speed 
of relative motion during the final approach. These constraints do not apply 
to the general case under consideration here, which more properly might be 
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thought of as a "phantom rendezvous" in which the target spacecraft is 
purely imaginary and therefore incapable of inflicting collision damage on 
the maneuvering vehicle. 

In relation to the targeting problem discussed in Section 3.2, adding 
an equality constraint on the phasing constant (t p ) of the target orbit has 
the effect of removing one degree of freedom in the problem solution and 
drastically reducing the number of possible transfer trajectories. This 
perhaps can best be illustrated by referring to Figure 10 in Section 3.2.4. 

If r Q is the -position of the initial impulse, then the maneuvering vehicle 
reaches the initial impulse point at a specific time t which is determined^ 
by the phasing constant of the initial orbit. Likewise, the phasing constant 

4 - 

of the target orbit determines the time t p at which the vehicle must reach 
the final impulse point defined by r p . The flight time in the transfer arc 
therefore must equal the difference in t p and t . If we restrict the dis- 
cussion for now to partial-revolution transfer arcs (i.e., where the transfer 
angle e is smaller than 360°), then the flight time varies monotonically 
between zero and infinity along the hyperbolic locus of allowable transfer 
orbit velocity vector coordinates. It follows that only one discrete point 
on the hyperbola — hence only one particular transfer trajectory — will 
satisfy the target orbit phasing constraint. In general then, for arbitrarily 
fixed impulse locations such as illustrated in Figure 10, it is not possible 
to satisfy a target orbit phasing constraint simultaneously with an equality 
constraint on the velocity increment magnitude of either {much less both) of 
the impulses. 

There is no expl icit solution for calculating the discrete point on the 
hyperbolic locus that will satisfy the flight time constraint discussed in the 
previous paragraph. There are a number of efficient iterative algorithms 
(e.g.. Reference 7) for calculating a fixed-flight- time transfer trajectory 


The term "spatial rendezvous" has sometimes been used to distinguish this 
problem from that of true rendezvous with a physical object; however, the 
terminology used above is thought to be less ambiguous. 

^Subject to discrete variations which are equal to integer multiples of the 
appropriate orbit period. 
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from one arbitrary point to another in an inverse-square gravity field. 

Most of these are based on Lambert's theorem regarding the relationship 
between flight time, the semimajor axis of the transfer orbit, and the 
lengths of the sides of the triangle illustrated in Figure 10. For this 
reason, this particular problem type (where the flight time and the impulse 
point locations are fixed) is generally characterized as “Lambert's problem" 
and the solution to such a problem is referred to as a "Lambert solution." 

When the initial and final orbits are fully defined (i.e., when p, e, 
i, Q , w, and t p are known for both orbits) as they are in this case, then 
the transfer trajectory flight time and -the locations of both impulses can 
be fixed unambiguously by assigning values to the impulse times t and t p . 
Since it is not possible to satisfy an equality constraint on aV at either 
impulse point when t and t p are assigned arbitrary values, we are reduced 
to the stratagem of systematically varying (i.e., scanning) the impulse 
times and computing a Lambert solution for each discrete combination of 
values in the hope that some particular combinations will produce the 
desired velocity increment magnitude(s) . 

The Two-Maneuver Spatial Rendezvous Processor (SR2MAN) of the prototype 
Mission Design and Analysis Subsystem (MDAS), which is described in Reference 
4, is configured to facilitate scanning operations of the type described in 
the preceding paragraph. In SR2MAN, the scan-control variables consist of 
one impulse time and a transfer-arc flight time. The specified impulse time 
may apply either to the first or the second maneuver, and the flight time 
may be either positive or negative. Based on the sign of the flight time, 
SR2MAN deduces whether the specified impulse time applies to the first or 
second maneuver, and computes the time of the other impulse accordingly. 

During a scanning operation, the SR2MAN processor is executed repeatedly 
under the control of the MDAS Monitor Program, according to user-supplied 
specifications that define the scope and resolution of the scan. The value 
of the impulse time and the flight time is fixed for each SR2MAN execution. 
After each execution, the Monitor program stores the SR2MAN output data in 
an array called a "data box" for subsequent processing, and increments the 
impulse time and/or the flight time before returning control to SR2MAN so 
that it can compute a new solution. 
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Each time it is executed, SR2MAN calculates a value for each of 27 
different scalar dependent variables that describe the nature of the trans- 
fer trajectory and of the associated maneuvers. After the completion of a 
scan, all of the output values that were generated during the scan reside 
in the data box mentioned in the preceding paragraph. Employing the Data 
Box Display Processor (DBDISP), the MDAS user can create two-dimensional 
quasi -graphical displays which summarize the behavior of selected dependent 
variables over the complete scan space. Figure 12 is an example of such a 
display, where the dependent variables are the magnitudes of the initial and 
the final velocity increment. Crude contours have been sketched in to illus- 
trate the approximate loci of solutions as they might appear for a particular 
combination of SRM stage AV capabilities. The solution of a singly-constrained 
or a doubly-constrained problem can be located as accurately as may be nec- 
essary (if it exists) by "zooming in" on the appropriate region of the display 
space in subsequent scans of lesser scope and greater resolution. 

If the flight time between the initial and the final impulse is suffi- 
ciently great, it may be possible to find multiple solutions to Lambert's 
problem by allowing the spacecraft to coast through one or more complete rev- 
olutions in the transfer orbit. Each complete revolution that is allowed 
increases the number of possible solutions by two, one of which is charac- 
terized by a larger transfer orbit semimajor axis than the other. An integer 
input variable IXPSOL is used in SR2MAN to determine which solution is desired. 
A value of 1 designates that the basic partial-rev solution is desired. 

To obtain a multi-rev solution the user sets IXPSOL = 2 M + L, where M is 
the number of complete revolutions to be allowed in the transfer orbit and 
where L has the value 0 or 1 depending on whether the transfer orbit 
with the larger or the smaller semimajor axis is to be selected. Of course, 
the value of IXPSOL is always fixed throughout any particular scanning 
operation. 

Although the SR2MAN processor as it is presently configured affords 
the basic capability (in conjunction with the DBDISP processor) for solving 
the phantom rendezvous problem as it applies to SRM stages, some modifica- 
tions are needed to make it fully effective. Using impulse time and flight 
time as independent variables for scan control produces -twofold and manifold 
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ambiguities regarding the orientation of the transfer orbit angular momentum 
vector. These ambiguities are essentially the same as those described in 
Section 3.2.1 for the case where Aj and Ay are considered as scan control 
variables. The SR2MAN logic was formulated originally for use in the design 
of trajectories for liquid propellant engines, where generally the objective 
is to minimize the velocity increment magnitudes. Accordingly, within the 
SR2MAN logic the ambiguities were resolved arbitrarily in favor of smaller' 

A V values. An alternate pair of scan-control variables, such as impulse 
time and wedge angle of the first maneuver, should be provided so that the 
manifold ambiguity associated with nodal transfers can be resolved fully. 

An auxiliary input variable (such as the variable H defined in Table 1) 
should also be provided so that the resolution of the twofold ambiguity that 
arises in every solution can be controlled externally. 
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4. RECOMMENDATIONS 


The analytic solution for two-impulse orbital transfer without a 
phasing constraint (Section 3.2) should be given priority with regard to 
software development. This is the most complicated of the trajectory design 
problems considered in this report, and probably the one that will have to 
be solved most' frequently in connection with Shuttle payload operations in- 
volving solid rocket motors. Preferably the software should be designed, 
along the same lines as SR2MAN 1 , to operate within the MDAS framework. This 
will allow the existing MDAS Monitor and the DBDISP processor to be used 
for scan control and data display. 

The SRM targeting requirements for interplanetary injection from a 
Shuttle parking orbit need to be analyzed in detail. For a given injection 
point in a fully-defined parking orbit and a specified arrival time at the 
target planet, all three components of the injection velocity impulse are 
uniquely determined (subject to very slight-variations with varying values 
of target parameters at the destination planet). Satisfaction of an equality 
constraint on the magnitude of the velocity increment will therefore be pos- 
sible only at certain discrete points in the parking orbit. By assuming 
that v ro is fixed (which is reasonable for analytical purposes) it appears 
from a preliminary study of the appropriate equations that an explicit solu- 
tion for the impulse points might be possible. However, the complexity of 
such a solution (if it exists) appeared to be such that it would be more 
properly treated as the subject of a separate report. 

Finally, with regard to the algorithms discussed in Section 2.3, a 
decision regarding possible development of software for simulating SRM 
powered-flight attitude maneuvers might profitably be delayed until details 
of the IUS design become more firmly established. 
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APPENDIX: 


PROBABLE IUS DESIGN CHARACTERISTICS 



The following excerpts from Request for Proposal F04701-76-R-0057 (dated 
21 January 1976) and its attachments are presented for the purpose of 
identifying probable design and operational characteristics of the Inter! 
Upper Stage. 



A.l OPERATIONAL PERIOD 

"The IUS system segment extends the capability of the Space Shuttle, 
so that DOD and NASA spacecraft can be delivered beyond low-earth orbits 
during the 1980-1985 time period." (RFP Attachment 1, p. 1) 

A. 2 COMPATIBILITY 

"The IUS system is primarily to be compatible with the Space Shuttle. 
However, potential compatibility of the IUS with a Titan booster will be 
investigated during the validation phase." (RFP, p. 10) 

"The IUS vehicle and ASE shall (a) be capable of remaining in the 
payload bay for a minimum of three orbital periods, and (b) j a1so be compatible 
with abort, descent, and post landing operations." (RFP Attachment 2, p. 9) 

"The IUS vehicle with attached spacecraft should be capable of being 
deployed from the Orbiter-mounted ASE within 20 minutes, starting at anytime 
after insertion into the nominal Orbiter parking orbit($)." (RFP Attachment 2, 
p. 30) 

A. 3 OVERALL DESIGN 

“The IUS system is to consist of expendable solid propellant stage and 
appropriate support subsystems." (RFP, p. 9) 

"The IUS vehicle comprises the expendable stages that leave the Orbiter 
and place spacecraft in either extended earth orbits or in planetary mission 
orbits, in identified configurations suitable for the missions. It includes 
the stage structure, the solid rocket motors, the reaction control system, 
the avionics for guidance, navigation and control, telemetry, tracking and 
command communications, instrumentation, data management equipment, and the 
required electrical power and electrical cabling. The IUS vehicle also in- 
cludes the stage to spacecraft interfaces, and the stage to airborne support 
equipment interfaces." (RFP Attachment 2, p. 7) 

"The conceptual design steady-state limit load factor for the IUS vehicle 
should be limited to -5.0 g's in the longitudinal axis for earth orbital mis- 
sions. The associated lateral load factors, other than vibrational loads, 
should be + 2.0 g's in the pitch and yaw directions. These total loads shall 
include thermally induced loads," (RFP Attachment 2, p. 38) 
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"Conceptual design limit load factors should be limited to -5 g's in 
the longitudinal axis for earth orbital missions, -10 g's for three-axis 
stabilized planetary spacecraft, and -15 g's for spinning spacecraft. The 
associated lateral load factors should be limited to +2.0 g's in the pitch 
and yaw directions. The total load should include thermal induced loads." 

(RFP Attachment 3, p. 8) 

"The IUS system shall provide, by a spin table, the capability to spin 
spacecraft up to 100 RPM." (RFP Attachment 3, p. 8) 

A. 4 MAIN PROPULSION 

"Propulsion designs shall identify motor family sizes including motor 
off-loadings, performance parameters, control authority and predicted system 
life." (RFP Attachment 1, p. 8) 

The IUS contractor will be required to determine "the minimum practicable 
stable of propulsive modules needed to insure operational flexibility in the 
implementation of the DOD, NASA, and non-NASA mission models." (RFP 
Attachment 1, p. 13) 

Movable nozzles and fixed nozzles with LITVC are to be considered for 
thrust vector control during SRM design trade studies. (RFP Attachment 6, 

P- 10) 

“The SRM total impulse shall vary by no more than + TBD per motor." 

(RFP Attachment 2, p. 49) 

A. 5 REACTION CONTROL SYSTEM 

" RCS operational functions shall include (a) velocity vernier for SRM 
impulse uncertainty, (b) velocity increments, as for spacecraft post separation 
maneuvers, (3) single or multi axis attitude maneuvers, (d) coarse or fine 
mode attitude control, and (e) powered flight roll control." (RFP 
Attachment 2, p. 49) 

Cold gas, monopropellant, and bi propellant systems are to be considered 
in RCS design trade studies. (RFP Attachment 6, p. 10) 
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A. 6 GUIDANCE, NAVIGATION, AND CONTROL SYSTEM 


"The IUS vehicle shall have stable attitude control, about the guidance 
commanded attitude, for all phases of coast and powered flight operations." 
(RFP Attachment 2, p. 27) 

"The IUS shall be capable of transmitting to the ground: (1) state 

vector and attitude data at final orbit injection of the 'IUS vehicle and 

spacecraft, and The IUS should also be capable of transmitting state 

vector and attitude data at transfer orbit injection." (RFP Attachment 2, 

P- 20) 

"The SGLS compatible radio links shall be capable of coherent turnaround 
of PRN modulation used to determine range and range rate at AFSCF ground 
stations." (RFP Attachment 2, p. 31) 

"The STDN compatible radio links shall be capable of coherent turnaround 
of the side tone ranging used to determine range and range rate at the STDN 
ground station." (RFP Attachment 3, p. 8) 

A. 7 PERFORMANCE 


The IUS vehicle will be capable of delivering a 5000 pound spacecraft 
to geosynchronous orbit or a 6000 pound spacecraft to a 12-hour orbit 
(h A = 21450 n mi, h p = 175/350/900 (min/nom/max) , i = 63.43°, u = 270°). 
(RFP Attachment 2, pp. 24 and 25) 

■ k 

The general position and velocity 3-sigma accuracies required at 
geosynchronous orbit injection are given below: 


Position (n mi) Velocity (fps) 


Tangential 

+ 66 

Normal 

+ 40 

Radial 

+ 50 

RSS 

+ 92 


+ 16 
+ 12 
+ 75 
+ 78 


(RFP Attachment 2, p. 26) 


These are designated as "Guidance and Navigation Accuracies" in the RFP 
attachment, but considering the context, it is believed that they are intended 
to represent the total allowable errors at the time of orbit injection. 
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"The IUS system shall be capable of placing single or multiple spacecraft 
in orbit by use of a single IUS vehicle in a single Shuttle launch. The IUS 
shall nominally place the spacecraft at a single orbital location (the space- 
craft being required to move to its final location). Relative phasing of the 
spacecraft shall not be required, except that For multiple deploy- 

ment, capability for a velocity increment of TBD shall be provided between 
each spacecraft deployment." (RFP Attachment 2, pp. 8 and 30) 

"The IUS vehicle shall be capable of performing maneuvers following 
separation which prevent direct IUS vehicle plume impingement on the space- 
craft and eliminate the possibility of subsequent re-contact with the 
spacecraft." (RFP Attachment 2, p. 31) 

For NASA geosynchronous missions, "The longitudinal placement for 
synchronous orbit spacecraft shall be +45°." (RFP Attachment 3, p. 6) 

"Ballasting, off-loading, pitch steering and other various energy 
management techniques should be considered." in motor sizing trade studies. 
(RFP Attachment 6, p. 10) 


*RFP Attachment 3 contains Design Reference Mission specifications for three 
specific planetary missions and one sun synchronous mission. 



